Theory of fractional vortex escape in a 0-k long Josephson junction 
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We consider a fractional Josephson vortex in an infinitely long 0-k Josephson junction. A uniform 
bias current applied to the junction exerts a Lorentz force acting on a vortex. When the bias 
current becomes equal to the critical (or depinning) current, the Lorentz force tears away an integer 
fluxon and the junction switches to the resistive state. In the presence of thermal and quantum 
fluctuations this escape process takes place with finite probability already at subcritical values of 
the bias current. We analyze the escape of a fractional vortex by mapping the Josephson phase 
dynamics to the dynamics of a single particle in a metastable potential and derive the effective 
parameters of this potential. This allows us to predict the behavior of the escape rate as a function 
of the topological charge of the vortex. 
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I. INTRODUCTION 

Josephson junctions with a phase drop of tt in the 
ground state (tt JJsi) are intensively investigated, as 
they promise important advantages for Josephson junc- 
tion based electronics^!^, and, in particular, for JJ based 
qubitsii^i^ii. Nowadays several technologies allow to 
manufacture such junctions: JJs with ferromagnetic 
barrie r^'^'^" , quantum dot JJ o^^i^^'^"^ and nonequilibrium 
superconductor-normal metal-superconductor JJsliii^. 

One can also fabricate O-tt long Josephson junctions 



(O-TT LJJs) 



16.17.18.19.20 



i.e., LJJs some parts of which be- 



have as junctions and other parts as tt junctions. The 
ground state phase fJ.{x) in such junctions will have the 
value of deep inside the 0-region, the value of tt deep 
inside the tt region and will continuously change from 
to TT in the Aj-vicinity of a O-tt boundary, where A,/ 
is the Josephson penetration depth. Such a bending of 
the phase results in the appearance of the magnetic field 
oc dfi/dx localized in the Aj-vicinity of a O-tt bound- 
ary and the supcrcurrcnts ±sin[^(a;)] circulating around 
it, i.e., one deals with a Josephson vortex. The total 
magnetic flux localized at the O-tt boundary is equal to 
$0/2, where $0 « 2.07 x 10~^^ Wb is the magnetic flux 
quantum. Therefore, such a Josephson vortex is called a 
semifluxon2ii22ii^. If the Josephson phase fi{x) deep in- 
side the TT region is equal to — tt instead of tt, one will 
have a localized magnetic flux equal to — <f>o/2 and a 
supercurrcnt of the vortex circulating counterclockwise 
(antiscmifluxon). Both semifluxons and antisemifluxons 
were observed experimentallySl and have been under ex- 
tensive experimental and theoretical investigation during 
the last decade . 

It turns out that instead of a 7r-discontinuity of the 
Josephson phase at O-tt boundary one can artificially cre- 
ate any arbitrary K-discontinuity of the phase at any 



point of the LJJ, and the value of k can be tuned 
clcctronically2i. As a result, in the ground state two 
types of vortices with the topological charges —k and 
— K -I- 27r (we assume that < k < 2tt) can be formed^. 
Such vortices are generalizations of semifluxons and an- 
tisemifluxons discussed above. They are stable only if 
their topological charge is smaller than or equal to 27r by 
absolute value. 

When one has a k vortex and applies a spatially uni- 
form bias current through the LJJ, the bias current exerts 
a Lorentz force, which pushes the vortex along the junc- 
tion. The direction of the force depends on the mutual 
polarity of the vortex and the bias current. The vortex 
exists just because it should compensate the phase dis- 
continuity, and. therefore, it is pinned in the vicinity of 
the discontinuity, i.e., it may bend under the action of 
the Lorentz force, but does not move away. Nevertheless, 
when the Lorentz force becomes strong enough, it tears 
off a whole integer fluxon out of a k vortex. The fluxon 
moves away along the junction, while a k — 27r vortex 
is left at the discontinuity. Further dynamics leads to 
the switching of the O-tt LJJ into the voltage state. This 
process was described for the flrst time for the case of a 
scmifluxon k = tt^ It takes place when the normalized 



bias current 7 
current o^i^ 



///co reaches the critical (depinning) 



7c (k) 



sin(K/2) 



k/2 



(1) 



where I^q = jcwL is the "intrinsic" critical current, which 
corresponds to the measurable critical current if k = 0; 
w is the LJJ width and L is its length. 

In the presence of quantum or thermal fluctuations, the 
escape process described above will take place with finite 
probability already at 7 < 7c. In this paper, we study 
the escape process of an arbitrary Josephson k vortex by 
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mapping the Joscphson phase dynamics to the dynam- 
ics of a single particle in an effective one-dimensional 
metastable potential. This allows us to predict the es- 
cape rates as a function of the bias current and of the 
vortex topological charge in the thermal and quantum 
domains. Our results can be directly compared with the 
experimental data that arc being obtained currently. 

II. MODEL 



The boundary conditions for fi{x,t) are 

^ix{-CO,t) = ^la;i+CO,t) =0. (6) 

At X' = the Josephson phase fi{x, t) and its derivative 
fix{x,t) are continuous, i.e., 

Ai(0+,t) = /i(0-,i) =Ai(0,i), (7a) 
/x^(0+,t) =/x,(0-,i) =//,(0,i). (7b) 



For our calculations we use dimensionless quantities. 
Lengths are measured in units of the Josephson length 
Aj, times are measured in units of where ujp is the 
plasma frequency, energies are measured in units of EjXj 
where Ej is the Josephson energy per length, and cur- 
rents are measured in units of the "intrinsic" critical cur- 
rent Ico of the Josephson junction. 

The dynamics of a fractional vortex in an infinitely long 
0-K Josephson junction with an applied bias current 7 is 
then described by the sine-Gordon equation (dissipation 
is neglected) 

^J'xxix, t) - ntt{x, t) - sin [ii{x, t) + kH{x)] = -7 , (2) 



where H{x) is given by 



Hix) 



for a; < 

1 for a; > 0. 



(3) 



This equation can be derived from the Lagrangian den- 
sity 

where the potential energy density U {fi, x) is given by 
U [fi,x] = 1 — cos + kH{x)] — 7^. (5) 



III. PROPERTIES OF THE STATIONARY 
SOLUTION 

The stationary solution ^0 (x) of Eq. ^ follows from 

fi'^ix) - U' [fio{x),x] = 

H'^{x) - sm[^io{x) + KH{x)]+-f = 0. (8) 
Without solving this equation we can derive some prop- 
erties of fJ,o{x) which will be useful for our calculations. 
Multiplying Eq. ([8]) with fiQ{x) in the regions a; < and 
a; > leads to 

1 2 

- 2 il^'oi^)^ + U [Mo(a;), x] = const. (9) 

Using the boundary conditions ^ for the stationary so- 
lution and the abbreviations 



^o(-oo); 
^o(+oo) + 



(10a) 
(10b) 
(10c) 



we obtain 



U [fio{x),x] - U [^o(-oo),-oo] cos(p_ - cos [ijlo{x)\ + 7 [(^_ - Mo(a:)] 

U [^j.q{x),x\ — U [^o(+oo), +00] = cosi^+ - cos [/io(a:) + k] + 7 [ip+ — Mo(a;) — . (H) 

These two equations allow us to express A*o(a^) in terms of fiQ^x) and ip± [or /io(±oo)]. Furthermore, since ij.q{x) and 
/Xo(x) are continuous at a; = 0, we have 

U [/io(0), 0-] - U ifioi-^), -cx)] = U [mo(0), 0+] - U [^ioi+^), +(X)] (12) 

which we can rewrite in the form 

cos (p+ — COS (p^ + J ((/?+ — (fi^ — k) — cos [ipo + k/2] — cos [(po — k/2] — —2 sin — sin ipo . (13) 

For X — > ±c» the phase /io(a;) approaches a constant value, which minimizes (stable stationary solution) or maxi- 
mizes (unstable stationary solution) the potential energy density U [/io(x), x], Eq. For stable stationary solutions 



a: < : - [fioix)f = 
a; > : i WQ{x)f = 
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wc find 

sin<^±=7, cos(/3± > =^ tp± = arcsin^ + 2n±TT , cos f± = ^^1 — 7^ , (14) 
whereas for unstable stationary solutions we find 

smip±=j, cos(p± < ip± ~ TT — avcsm"/ + 2n±7r , cos = — yl — 7^ . (15) 



Using tfiese results, Eq. ([T5)) can be reduced to 

7 (2n7r — k) = —2 sin — sin (po , (16) 

where n = 71+ — n_ is the number of fluxons already 
present in the system. Obviously, this condition cannot 
be fulfilled for arbitrary values of 7. At a critical value 7c 
we will find siniySp = ±1. For larger values of 7, Eq. (jl6p 
has no solution for ipo^ and the stationary solution fJ-oix) 
cannot exist. For n = 0, i.e., ip-^- = ip-, we obtain the 
critical current given in Eq. ([1]). 

IV. EIGENMODES 

The stability of the stationary solution /io(a;) can be 
analyzed with the help of the eigenmodes of the sine- 
Gordon equation ([2]). To find these eigenmodes we insert 
the ansatz 

fi{x,t)=fioix)+4'{x)e~'^' (17) 

into the sine-Gordon equation ^ and linearize it, as- 
suming I "0(2:) I ^ 1- Since fiQ{x) solves the stationary 
sine-Gordon equation, we obtain the Schrodinger equa- 
tion 

-i^'[{x) + U" [/io(a;), x] ?A„(a;) = 
= —^P'jKx) + cos[fio{x) + kH{x)]'4j„{x) = 
= u;lMx) (18) 

for the eigenmodes ipn(x), where the index n enumerates 
the eigenmodes. In this Schrodinger equation, the poten- 
tial is determined by the stationary solution fJ-o{x). Since 
the boundary conditions for fi{x, t) are already taken into 
account by the stationary solution fio{x), the boundary 
conditions for the eigenmodes ipnix) read 

V';(-c^)=V:.(+c^) = 0. (19) 

At a; = the phase ii(x,t) and its derivative /j,xix,t) are 
continuous, see Eq. ([7]). Since the stationary solution 
^o{x) and its derivative /j'q{x) are continuous at x = 0, 
ipnix) and ip'ni^) have to be continuous at a; = as well, 
i.e., 

V'n(0+) =M0-) = V-nlO); (20a) 
V';(0+) =^'n{0-)= C(0). (20b) 



As long as all eigenvalues of the Schrodinger equa- 
tion p8)) are positive, the stationary solution ^o{x) is 
stable. As soon as one eigenvalue becomes negative, 
the stationary solution fJ.o{x) is unstable. Therefore, at 
the critical current 7 = 7c the lowest eigenvalue, denoted 
by ujq, becomes zero. 

We now briefly show that at 7 = 7c the lowest eigen- 
mode "00(2;) is the derivative of the stationary solution, 
i.e., "ipoix) = C^'q{x). By taking the derivative of the 
stationary sine-Gordon equation ([8]) in the regions .t < 
and a; > we find 

- ^ /io(a;) + cos [tio{x) + kH{x)] fi'o{x) = . (21) 

Therefore, we have found the formal solution ipo{x) = 
C^q(x) of the Schrodinger equation with an eigen- 
value = 0. 

The boundary conditions (jl9p and the matching con- 
ditions for the eigenmodes introduce additional con- 
ditions for the stationary solution fio{x): 

?/;o(±oo) = CyUo(±oo) = ^ A'o(±oo) = 0, (22a) 

V'^(0+) = C^,'^{0+) = i''oio-) = c^i'iio-) 

^H'^i0+)=fi'^{0-). (22b) 

Using the stationary sine-Gordon equation Q and the 
abbreviations defined in Eq. (jlO[) . we can rewrite these 
conditions in the form 

sin(/3_ = sini^_|_ ~ 7, (23a) 
sin((^o — '*/2) = sin((y9o + k/2) 

K 

^ sin — cos ifiQ = 0. (23b) 

The first condition is fulfilled for any value of 7 [see 
Eqs. (fT4|) and ^5]) ]. whereas the second condition is only 
true for 7 = 7c, where we have sin 1^90 = ±1 and cos(po = 
0, see the discussion after Eq. p6)) . Therefore, tpo{x) = 
Cij,'q{x) is only an eigcnmode for 7 = 7c, see also Ref. l39l. 

V. EFFECTIVE POTENTIAL 

The phase fi{x, t) can be written in the form 

li{x,t)^ HQ (x) + ^ g„ (t) 0„ (x) . (24) 

n 

By inserting this expansion into the Lagrangian den- 
sity ^ and integrating over x we can derive a Lagrangian 
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for the mode amplitudes qn{t), which describes the mo- 
tion of a fictitious particle in many dimensions. For 
7 close to 7c the eigenfrequency ujq approaches zero, 
whereas the other eigenfrequencies remain finite. Around 
the minimum at qn = the potential will be "flat" in 
the direction of qq and "steep" along the other direc- 
tions. Therefore, we expect that at low energies, a par- 
ticle trapped in the minimum of the potential will move 
along go- Motivated by this simple picture, we only take 
into account the dynamics of the mode amplitude qo{t). 



To simplify the notation we denote the amplitude of the 
eigenmode ipo{x) by q{t). 

The Lagrangian for q(t) can be derived by inserting 
the ansatz 

^j,{x,t) fio{x) + q{t)iljQ{x) (25) 

into the Lagrangian density ([4]). Our simplified La- 
grangian reads 



-t-CX3 -t-OO +00 +00 

i = y Cdx=^q'^{t) J iIjq (x) dx - ^ J [hq{x) + q{t) ipoix)]^ dx - J U [no{x) + q{t) tpoix) , x] dx , (26) 

— oo — oo — oo — oo 

where £ and U are defined in Eqs. (|4]) and ([5]). Since we want to describe the escape of a particle from a metastable 
potential we only take in to account terms up to third order in q(t) and use the approximation 

U [fJ.o{x) + q{t) tpo{x),x] w U [^lo{x), x] + q{t) U' [^lo{x), x] V'o(2^) 

+ i qHt) U" [/io(a^), x] i;l{x) + \ q\t) U'" K(x), x] V'o'(x) . (27) 
2 o 



After omitting a constant term we obtain 

L - ^<z'(i) / ipS{x)dx-q{t) / [fi'aix) ti'a{x) + U' [tio{x),x\Mx)} 



+ 00 +00 

- \<l\t) J {Woixf + U"[fio{x),x]i;l{x)}dx^^q\t) J U'" [fio{x),x]^l{x) dx . (28) 

— OO —oo 

Performing two partial integrations and taking into account the boundary conditions for fio{x) and ^o(^) leads to 

+ 00 +00 

L = ^f{t) J ^^{x)dx-q{t) J {- fi'^^{x) + U'[fio{x),x]}Mx)dx 

— OO —oo 

- Iq^it) J {~i^'^{x) + U"[Mx),x]Mx)}i'o{x)dx-^q''{t) J U'" [Mx),x]^l{x)dx . (29) 

— oo — oo 

The linear term vanishes since Ho{x) is the stationary solution of the stationary sine-Gordon equation In the 
quadratic term we can use —iPq{x) + U" [fj,o{x),x] ipoix) = ujq'iPo{x), see Eq. (fT8|) . Therefore, the Lagrangian can be 
written in the form 

L = i Mg2(i) _ 1m^2 ^2(^) _ ^Gq'it), (30) 

where M and G are given by 

+ OC +CXD +00 

M= iPl{x)dx, G= U'" [noix),x]i;^{x)dx = - sin[^io{x) + KH{x)]i;l{x) dx . (31) 



I 

This Lagrangian describes the motion of a fictitious particle of mass M in the effective potential 

VM - ^Mujlq^ + ^Gq\ (32) 
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which can be characterized by the frequency loq for small 
oscillations around the minimum ai q = and a barrier 
height 



AV 



3G2 



(33) 



These two parameters determine the escape rates in the 
classical and in the quantum regime, see Sec. IVIII Note, 
that the frequency for small oscillations around the min- 
imum and the eigenfrequency of the lowest eigenmode, 
calculated from Eq. (|18p . are the same. 

The calculation of the present section leads to the fol- 
lowing procedure to determine the frequency too and the 
barrier height A^: 

1. For a given bias current 7 < 7c we solve the sta- 
tionary sine-Gordon equation ([8]) numerically and 
find the stationary solution j-io{x). 

2. For this stationary solution ^j.o{x) we solve the 
Schrodinger equation (|18p numerically and find the 
eigenmode ipo (x) and the corresponding eigenvalue 

3. Using the eigenmode i(jo{x) we calculate M and G 
numerically to find the barrier height AV. 



VI. ANALYTICAL APPROXIMATIONS 



A. Approximations for M and G 



The approximations for AI and G for 7 close to 7c are 
straightforward. Since M and G remain finite at 7 = 7c 
we replace M and G by their values at 7 = 7c and use 
the not normalized eigenmode i^o^x) ~ fJ-'oix), i.e., 



+00 



M 



G 



[fioix)] dx, 



(34a) 



-t-00 



sin [^J,o{x) + kH{x)] [^'^{x)] dx. (34b) 



These two integrals can be calculated analytically to a 
large extent. As shown in Appendix |3 the expressions 
for M and G can be written in the form 



</5o-k/2 

M = ±a/2 J ^ COS — cos <p -I- 7c (<y5- — 93) dkp ± \f2 J ^ cos — cos + 7c — V') '^'P (35) 



and 



G = 2 sin — sin (po cos 1,5+ -I- cos -I- 7c (<p+ -I- — 2(^o) 



(36) 



where and (^0 follow from the behavior of the to 

stationary solution for 7 = 7c at x = ±00 and x ~ 0, see 
Eq. pn]) . The upper sign for M applies to vortices with 
fi'o{x) > 0, whereas the lower sign applies to vortices with M 
fi'oix) < 0. 

For our examples in Sec. lVIIII we use vortices with q — 4sin- 



(7r+|K|)/2 

V2 / ^/ COS (fie — COS (fi + "fc i'Pc — (p)d(p,{38a) 



K 



COS tfc + 7c (<y3c - 7r/2) 



(38b) 



fiQ{— 00) = arcsin7, (37a) 
/^o(+oo) = arcsin7— K (37b) 

to satisfy Eq. In the limit 7 ^ 7c we therefore 

have = if- ~ arcsin7c = ipc which for positive bias 
currents implies (pQ — 7r/2. Moreover, h'q{x) is positive 
for —2tt < k < and negative for < k < 2tt. In this 
ease the expressions ([55]) and ([55)1 for M and G reduce 



B. Approximations for loq and AV 

Since ujq vanishes for 7 ^ 7c we cannot use its value at 
7c. Instead we use the stationary solution ^io{x) and the 
eigenmode ^^0(2;) = ^J'[){x) at 7 = 7c to evaluate the La- 
grangian ((29)) for 7 < 7c. In this case, the quadratic term 
in Eq. vanishes since U"[fio{x),x] does not depend 
on 7. The linear term, however, does not vanish since 
fJ-o{x) is not the stationary solution for 7 < 7c. Here we 
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have 



VII. ESCAPE RATES 



= -Mo (2;) + sin[/io(a;) + kH{x)] -^i = ^c- 7- (39) 

The remaining integral can easily be calculated and the 
approximate effective potential reads 



6 



where is given by 



+00 



+00 



(40) 



A^ ^ 

= /xo (+00) - /io (-00) = K. (41) 

For this potential wc find the frequency 



^o{x)dx= / y!Q{x)dx 



2GAfi 



Af2 



(7 - 7c) 



1/4 



(42) 



for small oscillations around the minimum and the barrier 
height 



2A/i 



(7 - 7c) 



-1 3/2 



3G^ 



(43) 



The first expression for AV" does not depend on M 
whereas the second expression for AT^'^'' agrees with 
Eq. ((33)l . except that now the approximate expressions 
for M, G, and loo are used. 



C. Point-like Josephson junction 

We would like to compare our results to the point-like 
JJ in the following sense. At given k, in experiments 
one sees a certain value of Ic and can calculate expected 
eigenfrequency, barrier height and escape rates using a 
short JJ model which ignores phase discontinuities and 
details related to the internal structure of the solution. 
We compare these results with our results for a fractional 
vortex. To obtain the results from the point-like JJ model 
we insert the values of 7 and 7c for a long Q-k JJ into the 
well known expressions for a point-like JJ^ and obtain 



Al/P* = 21 



7 

27c 



7 7 
— arccos — 



(44) 



and 



pt 



7c 



-1 1/4 



(45) 



where I = L/Xj is the normalized length of the Josephson 
junction. 



In Sec. |V] we have found that we can map the escape 
of a fractional vortex to the escape of a particle in a 
one-dimensional metastable potential which is character- 
ized by the frequency wq for small oscillations around the 
minimum and a barrier height A^. For such potentials 
approximate expressions for escape rates are available in 
the literature. As we use scaled quantities in the present 
paper we introduce the scaled temperature which mea- 
sures temperatures in units of EjXj/ks and the dimen- 
sionlcss parameter rj ~ hujp/ (EjXj) which plays the role 
of an effective h. 

In the classical regime, the escape of a particle from 
a metastable potential is due to thermal hopping. The 
escape rate (measured in units of ujp) for such processes 
is given by Kramers' formula^i^ 



- th 



-Av/e 

ZTT 



(46) 



where the prefactor p depends on the damping constant 
a of the system. For a/ujQ > 56/ {36AV) it reads 



P 




(47) 



If the damping constant a becomes too small, it has to 
be replaced by 



36aAV 
50LUO 



(48) 



In the quantum regime, the escape of a particle from 
a metastable potential is due to tunneling trough the 
energy barrier. In the semiclassical limit the decay rate 
of the ground state of a cubic metastable potential is 
given by^^ 



qm 



60 U)Q 



imAV f 36 AV\ 



(49) 



According to Ref. [43 for Ohmic damping the crossover 
between thermal hopping and quantum tunneling occurs 
at the temperature 



P 



T]UJo 

27r ' 



(50) 



where p is given by Eq. (|47p . 

Before we present our results we have a closer look at 
Eqs. (|46p and ()49p . in particular for small damping where 
we can use /? w 1. For both equations, the escape rates 
should be exponentially small, i.e., the exponents AV/9 
and 36AV/{5r]UJo) should be large. Using Eq. ([33]) we 
find that Kramers' formula, Eq. ((i^. can be apphed if 
luq satisfies 



'^o » 



2JVP 







^ 24 M3 a 



LOO > a/2. (51) 



The last two conditions allow us to use p « 1 . The semi- 
classical expression for the quantum mechanical escape 
rate (l49l) can be used for 



^0 » 



24 Af 3 



77 , Wo > 6'/??. 



(52) 



If the second condition is violated, the system will not be 
in the quantum mechanical ground state, and we have to 
take into account the quantum decay from exited states 
which is not included in Eq. For rjuiQ <C Al/ (many 

"bound" states in the mctastable potential) we may ther- 
mally average the decay rates from the ground state and 
the excited states. For details see Ref. |42|. 

The conditions introduced so far define lower bounds 
for iuo- On the other hand our simple model of a par- 
ticle moving in a onc-dimcnsional potential is only valid 
for 7 close to 7c where cuq becomes small. The main 
assumption we used to map the full problem to a one- 
dimensional problem was that lOq is much smaller than 
the other eigenvalues of the Schrodinger equation ^5)) . 
Therefore, we have to require Wq ^ w^, where ujf is the 
eigenvalue of the first excited state of the Schrodinger 
equation (fTS]). or the edge of the plasma band. 



VIII. RESULTS 

Our numerical results arc based on Eqs. ^ and ^TE\i . 
We solve these equations numerically for a symmetric 
junction with a length of 20Aj to emulate an infinitely 
long JJ and use vortices with /io(— 00) = /io(oo) + k = 
arcsin7 to satisfy Eq. (|14p . Two additional parameters 
are necessary to calculate escape rates: EjXj/ks is the 
conversion factor between the scaled temperature 9 and 
the temperature T, and r] = hujp/{EjXj) plays the role 
of an effective H. 

For our examples wc use the following JJ parame- 
ters: critical current density jc = lOOA/cm^, specific 
capacitance C = 4.2 /iF/cm^, junction width w = 1 ^m. 
This gives cjp = 27r • 42.8 GHz and EjXj = 78.4 meV = 
909 K. For these parameters, T and 9 are related via 
T = EjXj/kB -9 « 909 K • 0, and the value of 77 is 
77 = 2.3 • lO"'^. Furthermore, we assume that we are 
in the semiclassical limit, where we can apply Eq. 
and that we can use p w 1 in Eqs. ((46l) and (|50p . 
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FIG. 1: (Color online) The frequencies oj", ojq* and loq^ as a 
function of the bias current 7 for different values of k. The 
gray area indicates the plasma band, which was drawn by 
filling up the area above aji"(7). 



A. Comparison of different methods 



Before we present our results for escape rates, we first 
compare the approximations for luq and AV presented in 
Sec. IVII to the corresponding numerical values based on 
the single-mode approximation of Sec. |Vl For this pur- 
pose we calculate the eigenfrequency and energy barrier 
using three approaches: (a) single-mode approximation 
for 7 < 7c, Eq. ([TH]) for n = and Eq. denoted with 

superscript "nu"; (b) approximation at 7 = 7c; Eqs. (H^ 



and (|43|) . denoted with superscript "cr"; and (c) point- 
like J J formulas ([33]) and (|^5)) . denoted with superscript 
"pt" . 

The eigenfrequencies are shown in Fig. [TJ At k = the 
eigenfrequency " coincides with Wq * , while provides 
a very reasonable approximation. At larger values of k 
both and uj" provide good approximations to w™, 
but underestimate it a little. 



For our single-mode approximation to be valid, we 
have to make sure that the higher eigenfrequencies uin 
are much larger than loq. From Fig. [T]onc can see that 
this is not the case for small values of n where the eigen- 
frequency of a fractional vortex is close to the edge of 
the plasma band (shown in gray). Therefore, the single- 
mode approximation fails to describe the escape process 
in a long J J without discontinuities. On the other hand, 
Fig. [H shows the plasma band for an infinitely long JJ. 
For a JJ of finite normalized length / the plasma band 
consists of a set of discrete frequencies [j„, where the 
spacing between w„ is roughly inversely proportional to 
I. For moderate and especially for small values of I the 
difference between wo and the other eigenfrequency be- 
comes large, and the single-mode approximation works 
again even for k ^ (point-like J J formula) . 

The energy barrier calculated using different methods 
is shown in Fig. [21 For small values of k, AT^p* provides 
an excellent approximation to AF"" as expected in this 
limit, while ISV^^ quite overestimates the barrier. For 
large values of k (k = tt) the situation reverses: AT^'^'' 
approximates AT/"" very well, while AT/p' gives overes- 
timated values. The latter is expected as it is easier to 
activate a bent phase string starting the activation from 
some point than to move a fiat string simultaneously over 
the barrier. It seems that starting from k ^ 0.27r the 
values of /S.V'^'^ provide good approximations to AV™. 
Note, that we have chosen ^ = 20 to emulate an infinitely 
long JJ as the vortex solution is localized on the length 
scale 1 <C However, for /« the phase becomes flat 
and looses its localization, so that both AV"^ and AI^p* 
become a I, whereas AF" does not depend on /. In any 
case, for k — > the single (lowest) mode approximation 
does not work, so that one does not have to worry about 
the discrepancies in ISV in this limit. 



B. Quantum tunneling 

We use the numerically calculated values for the eigen- 
frequency Wg" and barrier height AV^"" and apply the 
semiclassical formula to these values to calculate 
quantum tunneling rates. The results are shown in Fig.[3l 
It turns out that the escape rate is larger, i.e., the vortex 
escapes easier, for larger values of k. This result is qual- 
itatively understandable because for k ^ the escape 
process reminds more and more the escape of the flat 
string from the metastable minimum and should scale 
with the JJ length. 



C. Thermal vs. quantum escape 

We use the numerically calculated values for the 
eigenfrequency ljq" and barrier height A]/"" and apply 
Eqs. with p = 1 and to calculate escape rates 
at different temperatures T, as shown in Fig. [Jfor k = tt. 
The dots represent the quantum escape rate. We remind 
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FIG. 2: (Color online) The energy barriers AV", A[/p* and 
AC/"" as a function of the bias current 7 for different values 
of K. 

that these rates were obtained using the Eqs. (|46l) and 
(Uni that are not valid very close to 7c (k), since in this 
case the escape rates are not exponentially small. Esti- 
mations using Eqs. ([STjl and ([52]) show that expressions 
(j46|l at r = 100 mK and expression (|49|1 become invalid 
for 7 > 0.9997c(7r). One can see that for T = 100 mK 
the two escape rates agree very well for 7 not very close 
to 1. We may therefore define a crossover temperature 
T* Si 100 mK — independent of 7 in the range of validity 
of (gni) and dH]). 
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FIG. 3: (Color online) Quantum escape rates Fqin as a func- 
tion of the bias current 7 for different values of k. 
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FIG. 4: (Color online) Thermal escape rates Fth as a func- 
tion of the bias current 7 for different temperatures. Dots 
represent the quantum escape rate Fqm- AH escape rates are 
calculated for k — n. 



pinning current in an infinitely long Josephson junction. 
By using a single (lowest) mode approximation, we have 
mapped the dynamics of an infinite dimensional system 
to the problem of a point-like particle escaping from a 
ID metastable cubic potential. For vanishing topologi- 
cal charge, the single mode approximation fails because 
the lowest eigcnmode is not well separated from the rest 
of the excitation spectrum. Thus, the lowest mode ap- 
proximation cannot be used to describe the escape in a 
conventional long J J (k = 0). 

In the region of validity of the single mode approxi- 
mation we have calculated the eigenfrequency and the 
barrier height numerically and analytically close to the 
depinning current. Then we have used the Kramers' 
formula and a semiclassical expression for thermal and 
quantum escape rates, respectively, to compare the es- 
cape rates of vortices with different topological charges 
and find the thermal-to-quantum crossover temperature. 
We have found that vortices with a larger topological 
charge escape easier. For typical experimental param- 
eters the crossover temperature lays in the range of 
100 mK as for many other JJ systems. These results can 
be directly compared to experiments that are in progress 
in the Tiibingcn group. 
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On the other hand, for the parameters used above 
Eq. dSni) with p = 1 and 7 ^ gives T* w 330 mK. 
According to Fig. [TJd) in the region of interest the eigen- 
frequency a;o(7r,7) is about 2 times smaller that a;(K, 0). APPENDIX A: INTEGRALS 
Thus Eq. dSni) predicts T* k. 160 mK in this region. 

In this appendix we evaluate the integrals for M and 
G. The calculations of this appendix are based on the 
unnormalized eigcnmode tpo{x) = fJ-Qix) which is only 
valid for 7 = 7c. Furthermore, we use the abbreviations 
(p- and ipo as defined in Eq. (|10p . 



IX. CONCLUSIONS 

We have investigated the thermal and quantum escape 
of an arbitrary fractional Josephson vortex close to its de- 
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With the help of Eq. ((TT|) the integral in the definition of G can be evaluated analytically for ^ = jc- 



+ CX3 



G 



sin [fJ-oix) + kH{x)] [^q{x)]^ dx 



sin [^J.oix)] [fJ-'oi^)]^ dx- I sin [^lo{x) + n] [^0(2^)]" d.^ 



m(o) 

sin 93 [cos 93_ — cos + 7c ~ <f)\ dip 
—2 J sin (cp + k) [cos (/J-). — cos {(p + k) + (<P+ ^ <y3 — k)] dv' 

Mo(0) 

fio — K./2 ip+ 

sin [cos — cos </? + 7c (</'- ~ v)] dip — 2 J sin [cos — cos (/? + 7c (</?+ ~ <y?)] dip 



cos^ — COS ip^ COS </5 — 7c sin ip — {'P- ~ cos 



2 

— cos"^ ip — cos cos ip — sin ip — 'jc {ip+ — ip) cos 



V- 

- V+ 

- ipo+K/2 



= — cos^(iy9o — k/2) — 27c sin((/3o — k/2) + 2 cos(</?o ^ k/2) [cos ip^ + 7c('^- — <,5o + k/2)] 
— cos^ + cos^ ip+ + 27c(sin cp„ — sin + cos^((/3o + k/2) — 27c sin(</Jo + '^/S) 
—2 cos((/?o + k/2) [cos + 7c(</'+ — ipa — '^/2)] • 
Using Eqs. (fT3|) and ([23|)we finally arrive at 



G = 2 sin — sin iy9o 



cos iy9+ + COS (y3_ + 7c + </3_ - 2ipo) 



In a similar way we can derive an expression for M. With the help of Eq. (|lip we obtain 



+00 



+00 



M 



[^I'oix)] dx^ Wo{x)] dx+ ifJ-'oix)] dx 



(Al) 



(A2) 



Mo(0) 



±V2 J Vc 



'cos(/7_ — cosiy9 + 7c (iy9_ — 1^9) dtp ± -\/2 y ■\/ COS ip+ — cos + k) + 7c {ip+ — ip ~ k) dip 

/Jo(-oo) /io(0) 



±V2 J V' 



cos — COS 1^9 + 7c (<rf5- — dip ± \f2 J \/ cos ip^ — cos <y9 + 7c (<p+ — ip^ dip . 

IP0 + K./2 



(A3) 



The upper sign applies to vortices with /io(x) > whereas t he lower sign applies to vortices with /Xn(x) < 0. 
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